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ABSTRACT 

A finite-volume Navier-Stokes code was used to simulate the shock-reflection and nozzle starting 
processes in an axisymmetric shock tube and a high Mach number nozzle. The simulated nozzle 
starting processes were found to match the classical quasi-one-dimensional theory and some features 
of the experimental measurements. The shock-reflection simulation illustrated a new mechanism 
for the driver-gas contamination of the stagnated test gas. 
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1 INTRODUCTION 

Free-piston driven reflected-shock tunnels have been one of the most successful types of impulse 
facilities [1] and have been able to provide a versatile range of hypervelocity test conditions with 
reasonable test-time durations. However, little is known of the detailed flow processes that are 
involved with their operation because, in the past, the unavailability of large computing resources 
has limited most analyses to either steady-state in two dimensions (as in nozzle design [2]) or to 
time- dependent simulation in one dimension (as in shock tube performance estimates [3] or nozzle 
starting estimates [4]). Recent studies have begun to investigate multidimensional transient effects 
such as the shock interactions in the throat region [5], axisymmetric nozzle starting times [6] and 
flow establishment over aerodynamic models [7]. This paper will study, via time-accurate numerical 
simulation, the shock reflection and the nozzle starting process in an axisymmetric shock tube and 
a high Mach number nozzle. 

1.1 REFLECTED-SHOCK TUNNEL OPERATION 

The principal features of a free-piston driven shock tunnel, along with an approximate wave diagram, 
are shown in Fig. 1. The driver (or compression) tube, which initially contains low pressure driver 
gas downstream of the piston, and the shock tube which contains the test gas, are separated by 
the (strong) primary diaphragm. Attached to the downstream end of the shock tube is the facility 
nozzle whose throat is significantly smaller than the inside diameter of the shock tube. The subsonic 
portion of the nozzle effectively closes the downstream end of the shock tube and forms the shock 
reflection region. The supersonic portion of the nozzle empties directly into a test section and dump 
tank which is evacuated to an initial pressure of approximately 30 Pa. The test gas is retained in 

the shock tube by a (thin mylar) secondary diaphragm. 

The first stage of operation consists of the launch of the piston and its acceleration along the 
compression tube. The driving force is supplied by compressed air from a reservoir. The driver gas 
in front of the piston is compressed and eventually bursts the primary diaphragm. After diaphragm 
rupture, the driver gas expands into the shock tube and shock-compresses the test gas before it. 
The primary shock wave travels the length of the shock tube, reflects from the closed end, and 
brings the test gas to rest in the nozzle supply region. 

Upon shock reflection, the light secondary diaphragm bursts and some of the test gas following 
ffie primary shock expands through the nozzle throat into the divergent part of the nozzle. From the 
point of view of the nozzle, the shock tube is now a reservoir of (nearly) stagnant, high-temperature, 
ligh-pressure test gas. The Mach 8 contoured nozzle discussed in this paper is large enough for 
;he starting time to significantly reduce the available test time. Thus, it was hoped that a direct 
limulation of the multidimensional flow processes would help identify the mechanism(s) delaying 
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the nozzle starting. See [8] for the classical model of the nozzle starting process. 

The useful test time is terminated by the arrival of driver gas at the nozzle exit. Generally, 
the mechanism used to model the premature driver-gas contamination of the test gas has been 
based on the interaction of the reflected shock with the boundary layer along the shock-tube wall 
[9]. However, in the simulations of the shock reflection process, an alternative mechanism is seen. 
Shortly after shock reflection, a vortex develops near the tube centreline and travels upstream with 
the reflected shock. This vortex enhances the mixing the fluid in the middle of the tube rather than 
along the tube walls. 


2 NUMERICAL SIMULATION 

The computations reported here were performed with a finite-volume upwind code based on the ful! 
Navier-Stokes equations. The code is described briefly here but further detail (including notation) 
is available in [10]. 


2.1 Governing Equations 

For an axisymmetric flow (with y as the radial coordinate), the finite-volume formulation of th< 
Navier-Stokes equations may be expressed as 

+ SF U yF - yF ’ } dy ~b Is {yG ~ ix = Q ' (1 

where < • > indicates a volume average, IT is the volume per radian and S is a contour in th< 
(a:, y)-plane around the volume. The £/, F and G vectors are 
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Treating the viscous contributions in the form yr avoids any difficulties with the geometry singularity 
at y = 0. The effective source term is 
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To augment these equations, we use the equation of state for a calorically perfect gas 

P = p{ 1 - l) e ? ( 7 ) 


Sutherland’s viscosity expression (to obtain y = y(T)), Stokes’ hypothesis A — 3 y and a constant 
3 randtl number. 


2.2 Numerical Implementation 

The flow domain in the (re, y)-plane is discretized as a structured mesh of quadrilateral cells with flow 
>roperties stored at the cell centres. At each time step, the inviscid-flux vectors (2) are evaluated by 
irst applying a MUSCL interpolation scheme to obtain “left” and “right” states at the midpoints 
if the cell interfaces. A locally one-dimensional (approximate) Riemann solver is then applied to 
■btain the interface flow properties during the time step. The spatial derivatives used in the viscous 
lux vectors (3)-(5) are obtained at the cell vertices by applying the divergence theorem. The source 
erm (6) is evaluated at the cell centres. The line integrals in (1) are then evaluated using the 
ciidpoint rule and the solution advanced in time with an Euler scheme or a predictor-corrector 

cheme. 

1.3 Flow Geometry and Initial Conditions 

’he flow domain used in these calculations included a 1.5m-long section of the shock tube, a 
ransition region into the throat and a Mach 8 nozzle. Fig. 2 shows the shock-reflection region and 
ozzle throat while the coordinates for key points are given in Table 1. The nozzle expansion consists 
f a conical section from a: = 0 to x = 0.168m followed by a contoured wall which straightens the 
ow (with the aim of producing a uniform and parallel test flow near the nozzle exit plane). The 
ontoured wall is specified as a cubic spline with the knots given in Table 2. 

The grids were generated as single- block structured grids. That is, a single pair of i,j indices 
-as used to address cells throughout the domain. This resulted in some mesh distortion in the 
ransition region between the shock tube and the nozzle throat. Cases run on a range of grid 


3 


resolutions, however, indicated that the mesh distortion had little effect upon the computed results. 
Starting at the entrance to the throat, cells were clustered toward the nozzle wall in order to resolve 
the boundary layer. No attempt was made to resolve the boundary layer in the shock tube or in 
the shock-reflection region. In the axial direction, cells were clustered toward the nozzle throat. 

Although the actual test gas used in the facility was air and the temperatures were high enough 
for chemical effects to be significant, the calculations presented here were run with a constant ratio 
of specific heats. A value of = 1.3 was used in order to approximate some of the high temperature 
effects. The gas constant was R — 287 J/kgjK and the Prandtl number was Pr = 0.72. Note that 
the driver gas was not included in these simulations. 

To approximate conditions discussed in [6], the initial state of the gas in the shock tube was set 
to 

Pi = 75.0 x 10 3 Pa, pi = 0.9679 kg/m 3 , m = 0 m/s, 7\ = 270 K, e x = 0.2583 x 10 6 J/kg, 
and a post-shock condition of 

P 2 = 5.812 x 10 6 Pa, p 2 = 6.764 kg/m 3 , u 2 = 2254 m/s, T 2 = 2994 K , e 2 = 2.864 x 10 6 J/kg, 

was applied to the inflow plane (x = —1.5m). This resulted in a primary shock Mach numbe: 
Mi = 8.29 and speed - 2630 m/s. At t = 0, the primary shock was located at x = —0.101m 
For ideal reflection, the stagnation conditions (between the reflected shock and the closed end o 
the tube) would have been 

Ps = 51.05 x 10 6 Pa, T s = 6330 K, p s = 28.10 kg/m 3 , 

with a reflected-shock speed of Wr = 714 m/s. The total enthalpy in the stagnation region wa 
= 7.87 x 10 e J/kg. If the gas was allowed to expand isentropically to sonic conditions, th> 
nominal throat conditions would have been 

P, = 27.86 x 10 6 Pa, p* = 17Mkg/m 3 , u* = 1433 m/s, T, = 5504 K, e, = 5.266 x 10 6 J/kg 

Bursting of the secondary diaphragm (located at x — 0) was modelled by running the simulatio 
in two stages. In the first stage, the primary shock was allowed to propagate to through the throa 
with pre-shock conditions applied as initial conditions throughout the nozzle expansion. As soon a 
the primary shock propagated into the nozzle expansion, the simulation was stopped and the flo 1 
state in the nozzle expansion (i.e. 0 < x < 2.187 m) was reset to 

Pi = 32.7 Pa, Pi = 3.833 x 10" 4 kg/m 3 , u t = 0 m/s, T, = 296 K, e< = 0.2832 x 10 6 J/kg . 

These conditions approximated the “evacuated” conditions found in the experimental facility. 
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3 RESULTS 

3.1 Nozzle Starting Process 

Fig. 3 shows the evolution of the density field on a 600 x 80 mesh. The simulation starts at t = 0 
with the primary shock propagating toward the throat. By t = 0.1 ms, the reflected shock is 
established and is bringing the test gas to rest upstream of the throat. Simultaneously, some of the 
stagnant test gas is expanding through the throat and into the divergent section of the nozzle. 

The essential features of the nozzle starting process are well described by quasi-one-dimensional 
model presented in [8]. As the primary shock travels down the nozzle it accelerates into the very low 
pressure gas already in the nozzle but, because of the diverging nozzle walls, it subsequently decel- 
erates. Test gas which accelerates through the nozzle throat following the primary shock, expands 
to a very high Mach number and is suddenly decelerated when it encounters the slower primary 
shock structure. An upstream-facing shock is thus formed and is swept downstream through the 
nozzle. Although the primary shock remains essentially planar, the contact surface between it and 
the upstream-facing shock undergoes a Rayleigh-Taylor (or Richmyer-Meshkov) type of instability 
and becomes quite distorted. The upstream-facing shock becomes slightly distorted but remains 

essentially coherent. 

Between the upstream-facing shock and the steady expansion being established near the throat 
there is an unsteady expansion, the upstream head of which is being swept downstream with velocity 
u-a. The unsteady expansion can be identified in Fig. 3 by the turning out of the contours between 
the upstream-facing shock and the nozzle throat. The unsteady expansion is the last wave to be 
expelled from the nozzle and by t= 1.2 ms the flow near the exit plane is approaching steady state. 

Fig. 4 shows the normalized Pitot profile (at late times) compared with experimental mea- 
surements from [6]. The key feature of the profiles is that they are all reasonably flat in the test 
core region. However, the computed boundary layer is too thin and consequently the average value 
of computed Pitot pressure in the core flow is lower than that measured in the experiment. The 
simulated and experimental Pitot pressure histories for three radial locations near the nozzle exit 
plane are shown in Fig. 5. The timimg of the passage of the upstream-facing shock and of the 
unsteady expansion closely agree but, the large disturbance seen in the experimental measurements 
for 0.7 ms < t < 1.0 ms at x ~ 0.006 is not present in the simulation results. 

3.2 Shock Reflection Process 

Although the scale of Fig. 3 is very small, it is clear that the reflected shock is not planar. At late 
times (i.e. t > 0.4 ms), this may be caused by inadequate grid resolution but, at earlier times, a 
complex structure is developed consistently for a range of grid resolutions. Fig. 6 shows a more 
detailed evolution of the density field in the shock-reflection region. To focus on events in the shock- 
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reflection (or nozzle-supply) region, the flow domain is restricted to —0.20m < x < 0.01m and is 
discretized as a mesh of 300 X 100 cells. The sequence begins at t = 0.030 ms with the reflected 
shock just beginning to move away from the flat end of the shock tube. Near the centreline, the 
primary shock continues to propagate downstream through the throat. At t = 0.040 ms, the 
reflected shock has moved further upstream and fluid which has been brought to rest behind the 
reflected shock begins to accelerate toward the tube axis and into the throat. Near the entrance to 
the throat (x ~ —0.025m, y ~ 0), a vortex is created and, in subsequent frames, continues to grow 
and follow the reflected shock. By t = 0.20 ms, the vortex has weakened and the reflected shock 
has become nearly planar. Note that this has not been observed in either experiment or simulation 
in two-dimensional flow. Fig. 7 shows the velocity field at an intermediate state and it can be seen 
that the vortex entrains fluid near the tube centreline well into the nozzle supply region. Hence, in 
the presence of a contact surface between driver gas and test gas (which is expected to be nearby 

for high-enthalpy operating conditions), this mechanism may cause premature contamination of the 
test gas. 

4 CONCLUDING REMARKS 


This paper has examined the shock-reflection and nozzle starting processes encountered in a reflected- 
shock tunnel. The simulation of the nozzle starting process did agree with the experimental data 
with respect to the timing of various events but did not identify the late-time disturbances seen in 
the experimental data. Obvious flaws in the simulation include, the lack of a turbulence model, 
and the omission of the driver gas and high-temperature chemical effects. 

The simulations did, however, identify a novel mechanism for the premature contamination 
of the test gas and may complement the traditional modelling of the interaction of the reflected 
shock with the shock-tube boundary layer [9]. This mechanism is based on the establishment and 
propagation of a vortex near the centre of the shock tube and just behind the reflected shock. 
Presently, this simulation has been done for only one operating condition, with a perfect gas, and 
with no boundary layers on the shock tube walls. More work needs to be done to investigate the 
mechanism in the presence of a driver-gas/test-gas interface and to find the range of operating 
conditions over which the mechanism is expected to be important. 
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label 

x, m 

r, m 

comment 

dl 

-1.5 

0.0 

inflow plane 


-1.5 

0.038 


b i 

-0.101 

0.0 

start grid distortion 

b 2 

-0.101 

0.038 



-0.063 

0.0 

turn grid 90° 

C 2 

-0.025 

0.038 


di 

-0.04262 

0.0 


d 2 

-0.025 

0.01762 

start radius into throat 

ei 

-0.015 

0.0 

start of constant area throat 

e 2 

-0.015 

0.00762 


ft 

0.0 

0.0 

end of throat 

f 2 

0.0 

0.00762 

start of conical expansion 


Table 1: Points defining the shock-reflection region and nozzle throat. 


H 


Vj , m 

comment 

i 

0.16800 

0.04951 

end of cone, dy/dx = 0,2493 

2 

0.45648 

0.10584 


3 

0.74488 

0.14208 


4 

1.03336 

0.16552 


5 

1.32176 

0.18104 


6 

1.61024 

0.19040 


7 

1.89872 

0.19520 


8 

2.18712 

0.19688 

end of nozzle wall 


Table 2: Knots for the cubic interpolating spline for the Mach 8 nozzle wall. 
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Figure 1: Reflected-shock tunnel operation, (a) Tunnel schematic; (b) x - t wave diagram. 



Figure 2: Geometry for the shock-reflection region. 
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Figure 4: Normalized Pitot pressure profile at x ~ 1.8m: O = experiment; o = simulation 
at t — 1.0ms; A = simulation at t = 1.2ms. 



t, ms 



t, ms 

Figure 5: Evolution of the normalized Pitot pressure at x ~ 1.8m: solid line = simulation; 
dashed line = experiment. Radial locations as shown. 
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Figure 6: Evolution of the density field (log 10 />) in the shock-reflection region. Times as 
shown. 
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